Jamming coverage in competitive random sequential adsorption of a binary mixture 



O 

o 

(N 
o 

CD 

Q 

m 



o 

CD 



CO 



I 

C 

O 

o 



> 

m 
o 
o 



o 



c 

o 
o 



X 



M. Kamrul Hassan ' , Jiirgen Schmidt 1 , Bcrnd Blasius 1 , and Jiirgen Kurths 1 
1 University of Potsdam, Department of Physics, Postfach 601553, D-14415 Potsdam, Germany 
University of Dhaka, Department of Physics, Theoretical Physics Division, Dhaka 1000, Bangladesh ^ 

(Dated: February 1, 2008) 

We propose a generalized car parking problem where cars of two different sizes are sequentially 
parked on a line with a given probability q. The free parameter q interpolates between the classical 
car parking problem of only one car size and the competitive random sequential adsorption (CRSA) 
of a binary mixture. We give an exact solution to the CRSA rate equations and find that the final 
coverage, the jamming limit, of the line is always larger for a binary mixture than for the uni-sized 
case. The analytical results are in good agreement with our direct numerical simulations of the 
problem. 
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The adsorption of particles onto a solid substrate is 
a common phenomenon and it is of wide interest in 
physics, chemistry, biology, and in other branches of sci- 
ence and technology. Examples include adsorption of gas 
molecules, colloidal particles, polymer chains, bacteria, 
proteins, and latex particles Jl|, |J ||. Due to its wide 
range of applications, adsorption has been studied exper- 
imentally, numerically, and analytically (see for example 
H for a review). 

The simplest model one can think of, still capturing 
the generic features of the adsorption phenomenon, is the 
kinetics of random sequential adsorption (RSA) of parti- 
cles of a fixed size on a d-dimensional substrate. The 
one-dimensional continuum version of the RSA process 
is popularly known as the random car parking (RCP) 
problem Q . In comparison to monolayer growth by RSA 
of monodisperse particles, very little attention has been 
given to the similar monolayer growth by two or more 
species of different size, despite the fact that the latter 
problem is much closer to the real life situation than the 
former. Nonetheless, there has been an increasing inter- 
est in the study of RSA of mixtures of different degree 
It is worth mentioning that the prob- 
lem is simplified in two extreme cases, namely, the RSA 
of uni-sized particles and that of mixture of particles of 
all sizes obeying a power-law size distribution JTcfl . In the 
former case the ultimate structure in the long time limit 
is described by the jamming coverage, whereas in the lat- 
ter case the resulting monolayer is uniquely quantified by 
the fractal dimension Df. The complexity arises in be- 
tween these two extreme cases when the mixture contains 
more than one species. The CRSA of a binary mixture 
has recently received new interest in the literature, where 
the properties of the jamming coverage obtained in ana- 
lytical models JllJ and in computer simulations f| have 
been discussed controversially. 

In this rapid communication, we consider the problem 
of the CRSA of a binary mixture on a one-dimensional 



substrate, i.e. the generalization of the classical RCP 
problem to two car species. To this end we develop an 
analytical model in terms of rate equations for the gap 
size distribution function and compute expressions for 
the jamming limit of both particle species, for a size ra- 
tio < 2. We compare our results to direct numerical 
simulations of the CRSA process. 

One time step of the CRSA process for a binary mix- 
ture is defined as: 

i) A particle is selected to be deposited on the sub- 
strate: With probability q the size of the particle 
is er and with probability p = 1 — q its size is ma, 
where m > 1. 

ii) Randomly a position is chosen on the substrate. 

Hi) The particle is deposited on the substrate if this 
position is not occupied by another particle and if 
there is no overlap with particles left and right to 
the position. 

iv) The process is repeated until no more gaps with 
size larger than a are left on the substrate. 

Our main results are as follows. The jamming coverage 
of the binary mixture always exceeds its monodisperse 
counterpart. The dynamics close to the jamming limit 
is governed by the dynamics of the smaller species. The 
total number density, as well as the number density of 
the smaller particles, gain their asymptotic values alge- 
braically, following Feder's law. The large particles reach 
their asymptotic coverage exponentially, with a decay 
constant ~ (m— l)q, multiplied by an algebraic prefactor 
t- 1 . 

We define the gap size distribution function 



P(x,t) = -N(x,t) 



(1) 



where N(x, t) is the total number of gaps with size x (say, 
in an interval x + dx) at time t on a substrate with length 
L. We scale length and time as 
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where a is the rate of particles brought to the substrate. 
In the limit L — > oo, keeping the particle inflow rate per 
unit length a/L fixed, the time evolution of the gap size 
distribution function is then described by the dimension- 
less rate equations 



dP(x, t) 
di 



dP(x,t) 
di 



dP{x,t) 
di 



-(x-{mp + q})P(x,t) (3) 

p OO pQO 

-2q P(y,t)dy + 2p P{y,t)dy, 

J x-\-l J x-\-m 



for m < x < oo 
-q(x-l)P(x,t) 



(4) 



2q P(y,t)dy + 2p P(y,t)dy, 

J x-\-l Jx+m 

for 1 < x < m 

2q P(y,t)dy + 2pl P(y,t)dy, (5) 

Jx+l 

for < x < 1 . 



These equations describe a RSA process with two par- 
ticle species competing for adsorption. The strength of 
competition is determined by q. The free parameters of 
the model are q and m. 

To solve the kinetic Eqs. (^]|), we seek a solution in 
the domain x > m of the form 



P(x,t) = A{t)e- [x - s) \ 



(6) 



where the average size of incoming particles is expressed 
as 

a = rap + q (7) 

and A(t) is to be determined. Substituting Eq. (^|) into 
Eq. (H) we obtain 



dlnAjt) 
Jt 



2q- 



t 



2p- 



t 



(8) 



At the beginning of the RSA process the number per unit 
length of gaps of any length x is zero. Also, initially there 
is only one gap of the size of the substrate. Solving Eq. 
(0) subject to these initial conditions, i.e. 



J dxP(x, 0) = 0, Hm o J dx xP 



gives 



where 



A(t) = t 2 F(t) 



(x,t) = l (9) 



(10) 



F(t) — e~ 2 ti du{q(l-e ")+p(l-e mu )}/u_ 

Thus, the solution for the gap size distribution function 
for gaps x > m is 



P(x,t) =t i F[t)e- { - x - 9 >, x > 



in 



(12) 



In the case m < 2 only the solution ( |l2| ) for x > m 
contributes to the integrals in Eq. ([|). Thus, in this case 
a solution in the domain 1 < x < m can be derived in 
the form 

P{x,t) = Bix,^-"^- 1 ^. (13) 
The result is (1 < x < m) 
P(x,t) = 2e-^ x ~ v > t 

dsF(s)s(qe imp - 1)s + p e - qms )e~ xps (lA) 



For m > 2 the yet unknown solution to Eq. con- 
tributes to the integrals in that equation as well, and 
thus, the integration of the equation is more difficult. 

From the knowledge of the gap size distribution func- 
tion in the domain x > 1 we find the total coverage Qt (t) 
by the adsorbing particles at time t 



(15) 



T (t) = 1 - / xP(x,t)dx. 



For < q < 1 the total coverage 9 T (t) = 6 s (t) + 9 L (t) 
has contributions by the small particles 0s(t) and large 
particles ^i(i). Practically, we first compute the time 
derivative 9r(t), using Eqs. (|3|-[|). Then we formulate 
the resulting expression in terms of P(x, t) in the do- 
mains x > m ( ^2] ) and 1 < x < m ( |l4| ) alone, by means 
of partial integrations. Thus, the solution to the rate 
equations for x < 1 is not needed. The contribution of 
the small particles 9s{t) to the coverage at time t is 



0s(t) 



F(s) 



2 ds- 

Jo qt+ps 



-s(p+mq) ^pg-f™- 1 ) 8 _|_ q^j 



;("»-!)«(—*) - s(qeS m -V s + p) + qt(e( m -^ s - 1)) 

+q [ e-( m -V q8 {l + (m-l)s}F(s)ds (16) 
Jo 

and the large particles contribute 

9 L (t)=pmf F(s)e- {m ^ 1 '> qs ds. (17) 
Jo 

We calculate the jamming limit of our model by numer- 
ical integration of Eqs. 
the classical RCP result 



j^ , |l7j ). For m = 1, we recover 



(t)= / F(s)ds, 



(18) 



which equals 9r — 0.748... (the Renyi-limit) for t — > oo. 

To test our analytic results we simulate the CRSA 
problem on a computer. Naturally, the simulations re- 
strict on finite substrate lengths L. However, for suffi- 
ciently large lengths L 3> a the effects of finite L are 
small. The algorithm we use follows literally the steps 
i) - iv) described in the introduction. The time scale in 



3 



the simulations is set by the number of particles brought 
to the substrate (successful or not). This discrete 'loop 
index' i relates to the time scale of the rate equations t 
as 



t 



L 



(19) 



Using this time scale and expressing all lengths in the 
simulation in terms of a, the simulational results can be 
directly compared to the solution to the rate equations 
and Eqs. (M, |l7|) for the jamming limit. 

In Fig. (pi) results for the jamming coverage are shown 
for the case m = 1.5, where the parameter q is varied. 
In general, we find an excellent agreement between sim- 
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FIG. 1: Coverage in the jamming limit vs q for m = 1.5, with 
the respective contributions from small and large particles. 
Symbols represent simulational results, lines the analytical 
solution ( [lq , |l7[ ). Simulations were carried out for a substrate 
length of 10000 small cars, each symbol is an average over 20 
independent realizations. 

ulation and theory. The total coverage is found to be 
always larger for binary mixtures (q ^ Q, q 1) than for 
the uni-sized case (q = or q = 1). Also, the contri- 
bution of the small cars to the jamming limit does not 
vanish as q — > 0. This is intuitively clear, since for a very 
small q first the large particles may reach their jamming 
limit Or = 0.748, while the small particles then have time 
to gradually fill the remaining gaps. In this sense the ad- 
sorption processes of large and small particles decouple in 
the limit q — > 0. The CRSA process is thus particularly 
effective in covering the substrate for small q and hence 
there is a sudden jump from 0x{q = 0) to 0r{q = + ). 

The highest coverages are obtained for small q and a 
large size ratio to of the binary mixture (see Fig. (Q)). In 
this case the large particles rapidly cover a fraction Or of 
the substrate, while the small particles then fill again the 
same fraction of the remaining gaps. Then the maximum 
possible coverage 9 max = 0r(2 — Or) = 0.937 may be 
reached in the limit m — >. The simulation with to = 20 
shown in Fig. (|J) is already very close to this limit. Our 
theoretical solution for the gap size distribution function 
restricts on the domain to < 2, thus no theory curve 
is plotted for the latter simulation. As q increases the 
coverage 0t decreases due to increasing strength of com- 
petition for adsorption between the two species. This 




FIG. 2: The total coverage vs q for various values of m. Di- 
amond symbols refer to simulations with a substrate length 
of 10000 small cars, and for the case m = 1.95 results from a 
simulation with L — 5000 are shown for comparison as square 
symbols. 



can be explained as follows. Every time a small particle 
wins in the competition it occupies less territory than a 
large particle would. Moreover, every occupation by the 
smaller specie creates an exclusion zone for further ad- 
sorption of larger particles. In the limit q — > 1 once 
again the two processes decouple but now the remaining 
gaps are smaller than the large particles and the symme- 
try is broken i.e. 0x{q — ► 0) ^ 0x(q — > 1). Note, that 
the CRSA model has no smooth transition in the limit 
to — > oo to the competitive adsorption model of a binary 
mixture of point particles and finite sized particles; the 
small particles of finite length will always occupy a frac- 
tion Or of the gaps between the large particles. Thus, 
the total coverage 0t is always larger than Or. Only if 
their length is indeed zero the small particles will not con- 
tribute to the coverage, and their only effect is to prevent 
the large particles from reaching their maximum cover- 
age Or. In this case the total coverage stems from the 
large particles alone (see Fig (|l|)), and it remains always 
below the Renyi- limit 9r fll4|. 
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FIG. 3: Time development of the coverage for q = 0.5, 
m — 1.5. Lines (solid, dashed, dotted) are theoretical results, 
symbols from simulations. 

The time development of the coverage is shown in Fig. 
(||) for the case q = 0.5, to = 1.5. Again a good agree- 
ment is found for theoretical and simulational results. 
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The contribution of the large particles to the total cover- 
age reaches its final value very rapidly. Asymptotically, 
Eq. @ yields 

/CO 
F(s)e-( m - 1)qs ds 

e -(m-l)qt 

~ , asi-^oo. (20) 

Thus, the asymptotic behavior of the coverage, reaching 
eventually the jamming limit, is dominated by the small 
particles. Indeed, the contribution of the small particles 
to the total coverage (Eq. (|l^)) approaches its final value 
algebraically as i _1 (see Fig. (Q)). Thus, we verified the 




TIME 

FIG. 4: The approach of 8s(t), and thus Jtr(t), to its asymp- 
totic value obeys Feder's law (equation (|2l|)). 

validity of Feder's law for this system 

e d (oo) - e d (t) ~ r^ d , (21) 

where d is the dimension of the substrate. 

The present work provides an exact theoretical basis 
on RSA of mixtures and removes the controversy between 
the numerical simulations presented by Meakin and Jul- 
lien for the two-dimensional CRSA of disks and that 
of Bonnier Jll| in Id. Meakin and Jullien report as well 
a larger coverage for the mixture. The dependence of our 
results on q and m is in complete qualitative agreement 
with the results reported in || ( see Fig. (3)). Further, 



Meakin and Jullien also confirm Feder's law for the small 
particle species, and also find an exponential approach of 
the coverage by the large ones. In the latter case, an 
algebraic pre- factor, as found here (Eq. (M)), might be 
hard to identify in their simulational data. 

In this context it is interesting to verify Palasti's con- 
jecture, stating that the jamming coverage 9d in d dimen- 
sions can be obtained from the one dimensional result, 
i.e. 

d (oo)«(0 1 (oo)) d . (22) 

Although the strict equality (as conjectured by Palasti 
|Q) is not valid, it is often found that the relation pro- 
vides a remarkably good estimate jl3| . Taking the values 
from figure 3 of the paper by Meakin and Jullien for the 
jamming coverages by a binary mixture of disks, we find 
a deviation of at most 5% from the estimates provided by 
our one-dimensional data and Eq. J2^ . The virtue of the 
present work is the excellent qualitative match of almost 
every espects of Id and 2d results. This confirms that Id 
model indeed is very useful in gaining information about 
heigher dimensions despite its simplicity. 

In conclusion we have shown that an analytical 
model for the competitive random sequential adsorption 
(CRSA) of a binary mixture in one dimension predicts 
a jamming limit that is always larger than for the uni- 
sized case. In the limits q = 0, q = 1, as well as m = 1 
we recover the classical result Or ; = 0.748 for the uni- 
sized RSA process. In all other cases, i.e. for true binary 
mixtures, we find a larger jamming limit in theory and 
simulation. In general, our analytical results are in ex- 
cellent agreement with direct numerical simulations. The 
time-asymptotic approach of the jamming limit is domi- 
nated by the contribution of the small particles, and we 
confirm the l/t behavior predicted by Feder's law in one 
dimension. The large particles reach their contribution 
to the jamming limit exponentially with an algebraic pre- 
factor. 
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